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By colliding two Bose-Einstein condensates we have observed strong bosonic stimulation of the 
elastic scattering process. When a weak input beam was applied as a seed, it was amplified by a factor 
of 20. This large gain atomic four- wave mixing resulted in the generation of two macroscopically 
occupied pair- correlated atomic beams. 
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Bose-Einstein condensation (BEC) in a dilute atomic 
gas produced a macroscopic population of the ground 
state wavefunction [||. Once BEC had been achieved, 
the initial well-defined quantum state can be transformed 
into other more complex states by manipulating it with 
magnetic and optical fields. This can result in a variety 
of time-dependent macroscopic wavefunctions Q , includ- 
ing oscillating condensates, multiple condensates moving 
relative to each other, an output coupler and rotating 
condensates with vortex lattices. Such macroscopically 
occupied quantum states represent classical matter-wave 
fields in the same way an optical laser beam is a classical 
electromagnetic wave. The next major step involves en- 
gineering non-classical states of atoms that feature quan- 
tum entanglement and correlations. These states are 
important for quantum information processing, sub-shot 
noise precision measurements [||, and tests of quantum 
non- locality. 

Quantum correlations in the BEC ground state have 
been observed in a BEC held in optical lattices |, |, |. 
The repulsive interactions between the atoms within each 
lattice site forces the occupation numbers to equalize, re- 
sulting in a number squeezed state. Alternatively, corre- 
lations in a BEC can be created in a dynamic or transient 
way through interatomic collisions. At the low densities 
typical of current experiments, binary collisions domi- 
nate, creating correlated pairs of atoms. Due to mo- 
mentum conservation, the pair-correlated atoms scatter 
into modes with opposite momenta in the center-of-mass 
frame, resulting in squeezing of the number difference 
between these modes [||, ^ 0. Our work is an imple- 
mentation of the suggestions in Refs. |||, ^. However, 
we use elastic scattering processes instead of spin fiip 
collisions to create pair correlations because the elastic 
collision rate is much higher than the spin fiip rate. This 
was essential to observe large amplification before further 
elastic collisions led to losses. 

Elastic scattering between two BEC's produces a col- 
lisional halo [^, where the number of atoms moving into 
opposing solid angles is the same, corresponding to num- 
ber squeezing. Once these modes are occupied, the scat- 
tering process is further enhanced by bosonic stimulation. 
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FIG. 1: Arrangement of laser beams to generate three atomic 
wavepackets: (a) Four Bragg beams intersected at the con- 
densate, (b) Two source waves and a small seed were created 
with Bragg beam pairs, (c) The four-wave mixing process 
amplified the seed and created a fourth wave. Both were sub- 
sequently read out with another Bragg pulse. The figures are 
projections on a plane perpendicular to the condensate axis. 
All wavepackets move within this plane. 



The onset of such an enhancement was observed in Ref. 
[||. In this paper, we report strong amplification, corre- 
sponding to a gain of at least 20. Based on a theoretical 
prediction [||, which drew an analogy to optical super- 
radiance we expected to obtain a highly anisotropic 
gain using our cigar-shaped condensate. However, this 
mechanism of mode selection proved to be irrelevant for 
our experiment, because atoms do not leave the conden- 
sate during amplification (see below). Instead, we pre- 
select a single pair-correlated mode by seeding it with 
a weak third matter wave, and observed that up to 40 
% of the atoms scattered into it. Because the scattered 
atoms are perfectly pair-correlated, the only fiuctuations 
in the number difference between the two beams stem 
from number fiuctuations in the initial seed. Therefore, 
an observed gain of 20 implies that we have improved 
upon the shot-noise limit by a factor of a/20, although 
this was not directly observed. Such a four-wave mix- 
ing process with matter waves had only been observed 
previously with a gain of 1.5 (lO[| 
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This experiment was performed with sodium conden- 
sates of ~30 million atoms in a cigar-shaped magnetic 
trap with radial and axial trap frequencies of 80 Hz and 
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FIG. 2: High-gain four wave mixing of matter waves. The 
wavepackets separated during 43 ms of ballistic expansion. 
The absorption images were taken along the axis of the 
condensate, (a) Only a 1% seed was present (barely visi- 
ble), (b) only two source waves were created and no seed, (c) 
two source waves and the seed underwent the four- wave mix- 
ing process where the seed wave and the fourth wave grew 
to a size comparable to the source waves. The gray circular 
background consists of spontaneously emitted atom pairs that 
were subsequently amplified to around 20 atoms per mode. 
The crosses mark the center position of the unperturbed con- 
densate. The field of view is 1.8 mm wide. 



20 Hz, respectively. Condensates had a mean field energy 
of 4.4 kHz, a speed of sound of 9 mm/s, and radial and 
axial Thomas-Fermi radii of 25 jam and 100 /im, respec- 
tively. The second condensate and the seed wave were 
generated by optical Bragg transitions to other momen- 
tum states. Fig. shows the geometry of the Bragg 
beams. Four laser beams were derived from the same 
laser that was 100 GHz red-detuned from the sodium D2 
line. The large detuning prevented optical superradiance 
[H. All beams propagated at approximately the same 
angle of ~ 0.5 rad with respect to the long axis of the 
condensate, and could be individually switched on and 
off to form beam pairs to excite two-photon Bragg tran- 
sitions [|l2| at different recoil momenta. 

The seed wave was created by a weak 20 /is Bragg pulse 
of beams with momenta pi and pa, which coupled 1-2 % 
of the atoms into the momentum state = p^ — p3, with 
a velocity ~15 mm/s. Subsequently, a 40 /is 7r/2 pulse 
of beams pi and p2 split the condensate into two strong 
source waves with momenta ki = and k2 = pi — P2 
(see Fig. |l|b), corresponding to a relative velocity of ~20 
mm/s. The four wave mixing process involving these 
three waves led to an exponential growth of the seed wave 
while a fourth conjugate wave at momentum k4 = ki + 
k2 — kg emerged and also grew exponentially (Fig. |l|c). 
The Bragg beams were arranged in such a way that the 
phase matching condition was fulfilled (the sum of the 
kinetic energies of the source waves (~ 11 kHz) matched 
the energy of the seed and the fourth wave). The effect 
of any energy mismatch on the process will be discussed 
later. The four- wave mixing process was analyzed by 
absorption imaging [Q. Fig. |^c shows the key result of 
this paper qualitatively: A small seed and its conjugate 
wave were amplified to a size where a significant fraction 
of the initial condensate atoms had been transferred into 
this pair-correlated mode. 
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FIG. 3: Absorption images after a readout pulse was applied 
to (a) the seed wave, and (b) the fourth wave. The thick 
arrows indicate the readout process. The readout pulse was 
kept short (40 /is), resulting in a large Fourier bandwidth and 
off-resonant coupling to other wave packets indicated by the 
narrow arrows. However, this did not affect the readout signal 
(atoms in the dashed box). 

To study this process, we applied "readout" beams p2 
and p4 for 40 /is, interrupting the amplification after a 
variable growth period between and 600 /is. (Turning 
off the trap after a variable amount of time is insufficient 
in this case because the density decreases on the time 
scale of the trapping period (1/(80 Hz)), while the am- 
plification occurs more rapidly.) The frequency difference 
between the two readout beams was selected such that 
a fixed fraction of the seed (fourth) wave was coupled 
out to a different momentum state k^ = kg + p4 — p2 
(k^ = k4 + p2 — P4) (Fig. |l|c). k^ or k^ did not experi- 
ence further amplification due to the constraint of energy 
conservation and therefore could be used to monitor the 
atom number in the seed (fourth) wave during the four- 
wave mixing process (Fig. ||). 

The growth of the 2 % seed and the fourth wave are 
shown in Fig. ||. As expected, the growth rates were 
found to increase with the mean field energy. Eventually, 
the amplification slowed down and stopped as the source 
waves were depleted. This is in contrast to Ref. [ p!Q| , 
where the mixing process was slow (due to much lower 
mean field energy), and the growth time was limited 
by the overlap time of the wavepackets. In our experi- 
ment, the overlap time was ^1.8 ms, whereas the growth 
stopped already after < 500 jas. 

A simple model describes the salient features of the 
process. The Hamiltonian of a weakly interacting Bose 
condensate is given by [0 

^ = 1^^<''-^—^ ««/.3«/.4 (1) 

where denotes the wavevectors of the plane wave states, 
m is the mass, V is the quantization volume, a^. is the 
annihilation operator and a = 2.75 nm is the scattering 
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length. If two momentum states ki and k2 are highly oc- 
cupied relative to all other states (with occupation num- 
bers Ni and A/'2), the initial depletion of ki and k2 can be 
neglected. Therefore, the only interactions are mean field 
interactions (self-interactions) and scattering involving ki 
and k2. In the Heisenberg picture, the difference between 



the occupations of the mode pairs An = a^^a,.^^ 



is time independent for any ki-\-K2 = ki+k2 . Therefore, 
the fluctuations in the number difference (An^) — (An)^ 
remain constant even though the occupations grow in 
time. The result is two-mode number squeezing. This is 
equivalent to a non-degenerate parametric arnplifier — the 
Hamiltonian for both systems are identical [Q. 

When calculating the occupations (al^^a^^), (^ils^'^s) 
and the correlation (^,^^a,^^)^ the relevant physical pa- 
rameters are 



Acj 



2m 2m 



mV 



N, (i = l,2) 



2m 



2m 



h 



h 



(2) 



where JI is the (geometric) average mean field energy of 
the two source waves, and hAco is the energy mismatch 
for the scattering of atoms from states ki and k2 to states 
hzi and hZ2' One obtains exponential growth for hzi and 
K:2 if p > hAco/A^ and the growth rate is given by: 
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For our initial conditions with s atoms in the seed wave 
kg and an empty fourth wave k4, the correlation (a\<id.\^^) 
start to grow as: 
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FIG. 4: Generation of pair-correlated atomic beams. The 
growth of (a) a 2 % seed and (b) it's conjugate fourth wave are 
shown for two different chemical potentials /i. The intensity 
of the waves are normalized to the intensity of the initial seed. 
The solid and dashed lines follow the initial growth according 
to Eq. (|) with growth rates of (170 /xs)"^ and (100 /is)"\ 
respectively. 



This leads to exponential growth of the occupation num- 
bers: 
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Eq. (|3|) and show that for a large mean field energy 
/I, Auj can be quite large without suppressing the four- 
wave mixing process. When Aco > Ajl/h^ one has to 
replace the hyperbolic sine functions in Eq. (^) with sine 

functions and r] in Eq. (||) with \J (Acc;/2)^ — {2'Jl/h)'^. 
The occupations in states k^ and k4 still grow initially, 
but then they begin to oscillate. The above solution also 
applies to initially empty modes with 5 = 0. 

We can estimate the maximum growth rate {2jl/h) for 
our experiment by using the average mean field energy 



across the condensate to obtain (53 /is)~^ for high and 
(110 /is)~^ for low mean field energy. The experimental 
data exhibit a somewhat slower growth rate of (100 /is)~^ 
and (170 /is)~^, respectively. This discrepancy is not 
surprising since our theoretical model does not take into 
account depletion and possible decoherence processes due 
to the finite size and inhomogeneity of a magnetically 
trapped condensate. We also observed that the angles 
of the Bragg beams and therefore the energy mismatch 
Acj, could be significantly varied without substantially 
affecting the four-wave mixing process, confirming the 
robustness (see Eq. (||)) of four- wave mixing. 

In addition to the four distinct wavepackets. Fig. |^ 
also shows a circular background of atoms that are scat- 
tered from the source waves ki and k2 into other pairs 
of initially empty modes hzi and ( + = ki + k2). 
The scattered atoms lie on a spherical shell in momen- 
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turn space centered at (ki +k2)/2 with a radius |k| close 
to |ki - ksl /2 and a width |Ak| - m^/^jf^JW/ |k|. As 
time progresses, the thickness |Ak| narrows due to the 
exponential gain. 

Eventually, the population of these background modes 
contribute to the depletion of the source waves. One 
can estimate the depletion time td of the source waves 
by comparing the total population in these modes to 
the original number of atoms. This sets a theoretical 
limit on the gain G = e'^^^^/^/4 given by 4G/Y^ln(4G) = 
\/2ii/ |k| a. For our geometry |k| a = 0.01 and G = 160. 
In our condensate of 3 x 10^ atoms, this maximum gain is 
achieved when all the atoms are scattered into the 9x10^ 
pair modes in the momentum shell. With a 1 % seed, the 
source waves are depleted earlier, leading to a maximum 
gain of 37, comparable to our measured gain of 20. 

In our experiment, we deliberately reduced the veloc- 
ity between the two source waves to twice the speed of 
sound in order to increase G and also the overlap time be- 
tween the two source waves. Under these circumstances, 
the thickness of the shell |Ak| becomes close to its ra- 
dius, accounting for the uniform background of scattered 
atoms rather than the thin s-wave halo observed in Ref. 
[||. For velocities around or below the speed of sound 
the condensate won't separate from the other waves in 
ballistic expansion. 

Once the amplified modes are populated, losses due to 
further collisions occur at a rate F ~ 87ra^n/z|k|/m per 
atom (n is the number density of atoms). In order to 
have net gain, the growth rate r] should be greater than 
F, which is the case since r]/T = l/|k|a = 100 ^ 1. Fur- 
thermore, we begin to lose squeezing when s + 1 atoms 
are lost from the mode pair that occurs approximately at 
a gain of e^^^/^/4 = l/|k|a. At this point, the condensate 
is already highly depleted. In our experiment however, 
the shell of amplified modes is so thick that it includes 
many of the modes into which atoms are scattered and 
increases the scattering rate by bosonic stimulation. Ide- 
ally, the atomic beams should separate after maximum 
gain is achieved. However, for our condensate size, the 
waves overlap for a much longer time and suffer colli- 
sional losses. This is visible in Fig. 2, where 40% of the 
atoms were transferred to the seeded mode pair, but only 
~ 10% survived the ballistic expansion. 

The collisional amplification process studied here bears 
similarities to the superradiant Rayleigh scattering of 
light from a Bose condensate where correlated 

photon-atom pairs are generated in the end-fire mode for 
the photons and the corresponding recoil mode for the 
atoms. However, there are significant differences between 
the two processes. In optical superradiance, the scat- 
tered photons leave the condensate very quickly, caus- 
ing only the recoiled atoms to maintain the coherence 
and undergo exponential growth. This physical situa- 
tion is refiected in the Markov approximation adopted 
in Refs. [^,0. In contrast, the atoms move slowly in 



collisional amplification, and the Markov approximation 
does not hold (although it was applied in Ref. The 
energy uncertainty AE = h/At for a process of dura- 
tion At gives a longitudinal momentum width of AE/v 
where v is the speed of light for photons or the velocity 
of the scattered atoms. This shows that optical super- 
radiance is much more momentum selective — the shell 
in momentum space is infinitesimally thin, and only the 
atomic modes with maximal overlap with this shell are 
selected. In contrast, the shell in collisional amplification 
is many modes thick and does not lead to strong mode 
selection. Moreover, in optical superradiance the light 
is coherently emitted by the entire condensate, whereas 
collisional amplification refiects only local properties of 
the condensate, because the atoms do not move signifi- 
cantly compared to the size of the condensate. Therefore, 
features like growth rate, maximum amplification, and 
even whether mode pairs stay squeezed do not depend 
on global parameters such as size or shape. 

In conclusion, we have observed high gain in atomic 
four-wave mixing and produced pair-correlated atomic 
beams. We have also identified some limitations for us- 
ing collisions to create such twin beams, including loss 
by subsequent collisions, and competition between other 
modes with similar gain. 
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